Time-odd mean fields in covariant density functional theory 
I. Non-rotating systems. 
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Time-odd mean fields (nuclear magnetism) are analyzed in the framework of covariant density 
functional theory (CDFT). It is shown that they always provide additional binding to the binding 
energies of odd-mass nuclei. This additional binding only weakly depends on the RMF parametriza- 
tion reflecting good localization of the properties of time-odd mean fields in CDFT. The underlying 
microscopic mechanism is discussed in detail. Time-odd mean fields affect odd-even mass differ- 
ences. However, our analysis suggests that the modifications of the strength of pairing correlations 
required to compensate for their effects are modest. In contrast, time-odd mean fields have profound 
effect on the properties of odd-proton nuclei in the vicinity of proton-drip line. Their presence can 
modify the half-lives of proton-emitters (by many orders of magnitude in light nuclei) and affect 
considerably the possibilities of their experimental observation. 
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I. INTRODUCTION 

The development of self-consistent many-body theo- 
ries aiming at the description of low-energy nuclear phe- 
nomena provides necessary theoretical tools for an ex- 
ploration of the nuclear chart into known and unknown 
regions. Theoretical methods (both relativistic and non- 
relativistic) formulated within the framework of density 
functional theory (DFT) and effective field theory (EFT) 
are the most promising tools for global investigation of 
the properties of atomic nuclei. The DFT and EFT con- 
cepts in nuclear structure models have been extensively 
discussed in a number of recent articles [2, H, E|. The 
power of the models based on these concepts is essen- 
tially unchallenged in medium and heavy mass nuclei 
where 'ab-initio' type few-body calculations are compu- 
tationally impossible and the applicability of spherical 
shell model is restricted to a few regions in the vicinity 
of the doubly shell closures. 

The self-consistent mean-field approach to nuclear 
structure represents an approximate implementation of 
Kohn-Sham density functional theory (DFT) 0, H, H, 0| , 
which is successfully employed in the treatment of the 
quantum many-body problem in atomic, molecular and 
condensed matter physics. The DFT enables a descrip- 
tion of the nuclear many-body problem in terms of a uni- 
versal energy density functional, and mean-field models 
approximate the exact energy functional, which includes 
all higher-order correlations, with powers and gradients 
of ground-state nucleon densities. Although it models 
the effective interaction between nucleons, a general den- 
sity functional is not necessarily related to any nucleon- 
nucleon (NN) potential. By employing global effective 
interactions, adjusted to reproduce empirical properties 
of symmetric and asymmetric nuclear matter, and bulk 
properties of some spherical nuclei, the current genera- 
tion of self-consistent mean-field methods has achieved a 
high level of accuracy in the description of ground states 
and properties of excited states in arbitrarily heavy nu- 



clei, exotic nuclei far from /3-stability, and in nuclear sys- 
tems at the nucleon drip-lines (see Refs. [1, [H, [l(| and 
references therein). 

The self-consistent methods (such as Hartree-Fock 
(HF) or Hartree-Fock-Bogoliubov (HFB)) based on zero 
range Skyrme forces or finite range Gogny forces are fre- 
quently used in nuclear structure calculations [j| [ill ]. 
These approaches represent non-relativistic energy den- 
sity functionals based on the Schrodinger equation for 
many-body nuclear problem [8[. 

On the other hand, one can formulate the class of rel- 
ativistic models based on the Dirac formalism, which 
can generally be defined as covariant density function- 
als (CDF) @. These models, such as quantum hadro- 
dynamics (QHD) 0, EH, are based on concepts of 
non-renormalizable effective relativistic field theories and 
DFT, and they provide a very interesting relativistic 
framework for studies of nuclear structure phenomena 
at and far from the valley of {3— stability Relativis- 
tic mean field (RMF) models Jl2| are analogs of the 
Kohn-Sham formalism of DFT [g], with local scalar and 
vector fields appearing in the role of local relativistic 
Kohn-Sham potentials [l], [l3j]. The exact energy density 
functional is approximated with powers and gradients of 
auxiliary meson fields or nucleon densities. The EFT 
building of the energy density functional allows error es- 
timates to be made, provides a power counting scheme 
which separates long- and short-distance dynamics and, 
therefore, removes model dependences from the self- 
consistent mean field approach jl4{ . In the description 
of nuclear ground states and properties of excited states 
the self-consistent mean-field implementations of quan- 
tum hadrodynamics, the relativistic Hartree-Bogoliubov 
model (RHB) and the relativistic (quasiparticle) random 
phase approximation (RQRPA) and their subversions, 
are employed 

The mean field is a basic concept of every DFT. One 
can specify time-even and time-odd mean fields (Tol . [l|| 
dependent on the response of these fields to the action 
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of time-reversal operator. The properties of time-even 
mean fields in nuclear density functionals are reasonably 
well understood and defined [1, This is due to the 
facts that (i) many physical observables such as bind- 
ing energies, radii etc. are sensitive only to these fields, 
and (ii) the model parameters are fitted to such physical 
observables. 

On the other hand, the properties of time-odd mean 
fields, which appear only in nuclear systems with bro- 
ken time-reversal symmetry, are still poorly understood. 
However, it is already known that these fields are impor- 
tant for proper description of rotating nuclei 15 L jUfT^]. 
band terminations Ell , magnetic moments [20||, large 
amplitude collective dynamics [2l[, fussion process [22 |. 
the strengths and energies of Gamow- Teller resonances 
[23j |. binding energies of odd-mass nuclei [HI, I2H I2II 
and the additivity of angular momentum alignments 1271 . 
They also may play a role in the N = Z nuclei [24l |28| 
and affect the definition of the strength of pairing corre- 
lations [HHH. 

There was a dedicated effort to better understand 
time-odd mean fields in the framework of the Skyrme en- 
ergy d ensity functional (EDF) theory (see Refs. [HI, [l7], 
ll8L |23| and references therein) . On the contrary, much 
less attention has been paid to these fields in covariant 
density functional theory (CDFT) [H, M, HE HI • The 
current manuscript aims at better and systematic under- 
standing of time-odd mean fields and their impact on 
physical observables in non-rotating nuclei in the frame- 
work of the RMF realization of the CDFT. The results 
of the study of these fields in rotating nuclei will be pre- 
sented in a forthcoming article [30( which represents a 
continuation of the current investigation. 

The manuscript is organized as follows. The cranked 
relativistic mean field theory and its details related to 
time-odd mean fields are discussed in Sect.fll] Section UTT1 
is devoted to the analysis of the impact of time-odd mean 
fields on binding energies of odd-mass nuclei. The mass 
and particle number dependences of this impact and their 
connections with odd-even mass scatterings are also con- 
sidered. The microscopic mechanism of additional bind- 
ing in odd-mass nuclei induced by time-odd mean fields 
is analyzed in Sect. IIV1 The impact of time-odd mean 
fields on the properties of proton-unstable nuclei is stud- 
ied in Sect. |Vl Section [VI] considers how time-odd mean 
fields modify the properties of odd-odd nuclei. Finally, 
Sect. IVIII contains the main conclusions of our work. 



II. THEORETICAL FORMALISM 

The results presented in the current manuscript have 
been obtained using Cranked Relativistic Mean Field 
(CRMF) theory [HI M, HI]. This theory has been suc- 
cessfully employed for the description of rotating nuclei 
(see Ref. [3J] and references therein) in which time-odd 
mean fields play an important role, but it is also able to 
describe the nuclear systems with broken time-reversal 



symmetry in intrinsic frame at no rotation. In this the- 
ory the pairing correlations are neglected which allows 
to better isolate the effects induced by time-odd mean 
fields. This is because time-odd effects also play a role 
in pairing channel [2l|. The CRMF computer code is 
formulated in the signature basis. As a result, the break- 
ing of Kramer's degeneracy of the single-particle states is 
taken into account in a fully self-consistent way. This is 
important for an accurate description of time-odd mean 
fields in fermionic channel (see Sect. IIV[) . Most impor- 
tant features of the CRMF formalism related to time-odd 
mean fields are outlined below (for more details see Refs. 
[3l], HH]) for the case of no rotation (rotational frequency 
n x = 0). 

In the Hartree approximation, the stationary Dirac 
equation for the nucleons in the intrinsic frame is given 
by 



(1) 



where hp is the Dirac Hamiltonian for the nucleon with 
mass m 

h D = a{-iV -V{r)) + V (r) + (3(m + S(r)). (2) 

It contains the average fields determined by the mesons, 
i.e. the attractive scalar field S(r) 



S{r) = g a <j{r), 



(3) 



and the repulsive time-like component of the vector field 
Vo(r) 

I — T3 

V {r) = g u u) (r) + g p T 3 p (r) + e — - — A (r), (4) 



A magnetic potential V(r) 

I-T3 

V(r) = g^{r) + g P T 3 p{r) + e——A(r), 



(5) 



originates from the space-like components of the vec- 
tor mesons. Note that in these equations, the four- 
vector components of the vector fields u>^, p^, and 
are separated into the time- like (luq, po and Aq) and 
space-like [u> = (uj x , u y , u z ) , p = (p x , p v , p z ), and A = 
(A x ,A y ,A z )] components. In the Dirac equation the 
magnetic potential has the structure of a magnetic field. 
Therefore the effect produced by it is called nuclear mag- 
netism (NM) [3l|. 

The corresponding meson fields and the electromag- 
netic potential are determined by the Klein-Gordon equa- 
tions 

{-A + ml} a(r) = -g a [p n s (r) + (f s (r)] 

-g2<J 2 {r)-g^\r), (6) 
{~A + mi}u; (r) = g u [p^(r) + pg(r)], (7) 
{-A + ml} «(r) = g^3 n (r)+j p (r)} (8) 

{-A + TOp} po(r) = g p [p-{r)-pl{r)], (9) 

{-A + ml) p(r) = g p [j n {r)-f{r)\, (10) 
-AA (r) = epP(r), —AA(r) = ef(r), (11) 



3 



with source terms involving the various nucleonic densi- 
ties and currents 



N,Z 



P n v ' P (.r) 
r , P{r) 



5^(V'i(r)) t ^i(r) 1 

N,Z 

£(^(r)) f ^(r), 

i=l 



(12) 
(13) 
(14) 



i=i 



where the labels n and p are used for neutrons and pro- 
tons, respectively. In the equations above, the sums run 
over the occupied positive-energy shell model states only 
(no-sea approximation) 35]. Note that the spatial com- 
ponents of the vector potential A(r) are neglected in the 
calculations since the coupling constant of the electro- 
magnetic interaction is small compared with the coupling 
constants of the meson fields. 

The magnetic potential V(r) in the Dirac equation as 
well as the currents j n ' p (r) in the Klein-Gordon equa- 
tions do not ap pea r in the RMF equations for time- 
reversal systems [12J . Similar to the nonrelativistic case, 
their presence leads to the appearance of time-odd mean 
fields. Thus, we will use the terms nuclear magnetism 
and time- odd mean fields interchangeably throughout 
this manuscript. The magnetic potential is the contribu- 
tion to the mean field that breaks time-reversal symmetry 
in the intrinsic frame and induces non- vanishing currents 
j n ' p (Eq. HI])) in the Klein-Gordon equations (Eqs. (0 
ITOIO . which are related to the space-like components of 
the vector mesons. In turn, the space-like components of 
the vector lj and p fields form the magnetic potential ^ 
in the Dirac equation. Note that the current j n ' p change 
the sign upon the action of time-reversal operator [361 ] - 
The currents j n ' p have relativistic meaning since they are 
the products of the large (F) and small (G) components 
of the Dirac spinors tpi( r )- 

The spatial components of the vector u> and p mesons 
lead to the interactions between possible currents. For 
the w-meson this interaction is attractive for all combi- 
nations (pp, nn and pn-currents), and for the p-meson it 
is attractive for pp and nn-currents but repulsive for pn- 
currents. Within mean field theory such currents occur 
only in the situations of broken time-reversal symmetry. 

Note that time-odd mean fields related to NM are de- 
fined through the Lorentz invariance and thus they 
do not require additional coupling constants: the cou- 
pling constants of time-even mean fields are used also for 
time-odd mean fields. 

The currents are isoscalar and isovector in nature for 
the u) and p mesons (Eqs. (JSJ QUI)), respectively. As a 
consequence, the contribution of the p-meson to mag- 
netic potential is marginal in the majority of the cases. 
Thus, time-odd mean fields in the RMF framework de- 
pend predominantly on the spatial components of the 
u> meson. Neglecting the contribution of the p meson, 



one can see than only two parameters, namely, the mass 
rriuj and coupling constant g u of the lu meson define the 
properties of time-odd mean fields (Eqs. ([5][S] and [TO])). 
Table Q] clearly indicates that these parameters are well 
localized in the parameter space for all modern non- linear 
parametrizations of the RMF Lagrangian. This suggests 
that the parameter dependence of the impact of time-odd 
mean fields on the physical observables should be quite 
weak. Indeed, the analysis of terminating states in Ref. 
fl9| showed that time-odd mean fields are defined with 
an accuracy around 15% in the RMF framework. 



TABLE I: The and g^ parameters of most frequently used 
non-linear parametrizations of the RMF Lagrangian. 



Parametrization 


muj 




NL1 (Ref. [35]) 


795.359 


13.2846 


NL3 (Ref. [37]) 


782.501 


12.8675 


NL3* (Ref. [38]) 


782.60 


12.8065 


NL-Z (Ref. [39]) 


780.00 


12.9086 


NL-Z2 (Ref. [40] 


780.00 


12.9084 


NL-RA1 (Ref. [41]) 


783.00 


12.9212 


NLSH (Ref. [42]) 


783.00 


12.9451 


TMA (Ref. [43]) 


781.95 


12. 842+3. 191A- 04 



The total energy of the system is given in Refs. [3l], 
[33j . In order to facilitate the discussion we split it into 
different terms as [66] 



E f „ 



E, 



part 



+ E B 



-E a - 

Ecoul , 



E r 



NL 



El - El 



(15) 



where E part and E cm represent the contributions from 
fermionic degrees of freedom, while other terms are re- 
lated to mesonic (bosonic) degrees of freedom. In Eq. 



(16) 



is the energy of the particles moving in the field created 
by the mesons (e^ is the energy of i-th particle and the 
sum runs over all occupied proton and neutron states) 



E„ 



1 



g a / d 3 r a(r)[ff s (r)+p^r)}, 



(17) 



is the linear contribution to the energy of isoscalar-scalar 
CT-field 



EaNL = ~ / '/ '; 



i ff2 <J 3 {r) + ^g 3 cr A (r) 



(18) 



is the non-linear contribution to the energy of isoscalar- 
scalar cr-field 



El = \g u I </ :! r^-u(r.» Ken + />;:<»•)] 



(19) 



is the energy of time-like component of isoscalar-vector 
w-field 



1 



K = o 5p d J rpo(r-)K(r)-^(r)] 



(20) 



is the energy of the time-like component of isovector- 
vector p-field 



e: 



i 



d 3 ru:(r)[f(r)+j n (r)}, (21) 



is the energy of the space- like component of the isoscalar- 
vector w-ficld 



F s - 1 



9 P / d s rp{r)\j n {r)-f{r)}, 



(22) 



is the energy of the space-like component of isovector- 
vector p-ficld 



Ecoui = ^ e / tfVA (rK(r), 



is the Coulomb energy 
3 



41A- 1/3 MeV, 



(23) 



(24) 



is the correction for the spurious center-of-mass motion 
approximated by its value in a non-relativistic harmonic 
oscillator potential. 

The CRMF equations are solved in the basis of 
an anisotropic three-dimensional harmonic oscillator in 
Cartesian coordinates characterized by the deformation 
parameters (5$ — 0.3 ((3q — 0.5 in the case of superde- 
formed states) and 7 = 0° as well as oscillator frequency 
hjjQ = AlA^ 1 / 3 MeV. The truncation of basis is per- 
formed in such a way that all states belonging to the 
shells up to fermionic Np=V2 and bosonic 7Vb=16 are 
taken into account; numerical analysis indicates that this 
truncation scheme provides sufficient numerical accu- 
racy. Single-particle orbitals are labeled by [Nn z A]^l st9n . 
[Nn z A]Q are the asymptotic quantum numbers (Nilsson 
quantum numbers) of the dominant component of the 
wave function. The superscripts sign to the orbital labels 
are used sometimes to indicate the sign of the signature 
r for that orbital (r = ±i). The majority of the calcu- 
lations are performed with the NL3 parametrization [37| 
of the RMF Lagrangian. 

In order to investigate the impact of NM (time-odd 
mean fields) on physical observables, the CRMF calcu- 
lations are performed with and without NM (the later 
will be further denoted as WNM). This is only the way 
in which the effects of time-odd mean fields can be stud- 
ied and as such it is frequently used in the DFT stud- 
ies, both in relativistic and non-relativistic frameworks 
[II [H GJ, [H H, [H [H OH ■ One should however keep 
in mind that if time-odd fields are neglected, the local 
Lorentz invariance (Galilean invariance in non-relativistic 
framework [l5|,[47|]) is violated. The inclusion of time-odd 
mean fields restores the Lorentz invariance. 
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FIG. 1: (Color online) The impact of NM on binding en- 
ergies of light odd-mass nuclei. The calculations have been 
performed with the NL3 parametrization of the RMF La- 
grangian. They cover the nuclei from the proton-drip line 
up to the neutron-drip line. 



The calculations of odd and odd-odd nuclei require 
that only one state among two time-reversal counter- 
part states is occupied in proton or/and neutron subsys- 
tem^). This is achieved in the CRMF code by occupying 
one signature (r = +i or r = —i) of the single-particle 
state, and keeping opposite signature empty. We will 
call the occupied signature as blocked state in order to 
simplify the discussion. 



III. BINDING ENERGIES IN ODD MASS 
NUCLEI 

The time-reversal invariance is conserved in the ground 
states of even-even nuclei. The nucleon states are then 
pairwise degenerated, and the contribution of the state 
to the currents cancels with the contribution of its time- 
reversed partner. Time-odd mean fields reveal them- 
selves in odd- and odd-odd mass nuclei and in two-(multi- 
)quasiparticle states of even-even nuclei. This is because 
unpaired (odd) nucleon breaks time-reversal invariance 
in intrinsic frame and produces the contribution to the 
currents and spin. In this case, the Kramer's degeneracy 
of time-reversal partner orbitals is also broken. 

The modifications of the binding energies and quasi- 
particle spectra are most important issues when consid- 
ering time-odd mean fields in non-rotating systems. The 
binding ener gies are important in nuclear astrophysics 
applications [481 ] . and their modifications due to time- 
odd mean fields may have considerable consequences for 
the r- and rp-process abundances. They are also one of 
the most valuable input data for the parameter fit of the 
nuclear structure models. Thus, it is important to un- 
derstand the influence of time-odd mean fields on bind- 
ing energies of odd- and odd-odd mass nuclei especially 
in the context of mass table fits |49l| . With current focus 
on spectroscopic quality DFT [501 ] . a knowledge on how 
the time-odd mean fields influence the relative energies 
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FIG. 2: (Color online) The same as in Fig. [T]but for the re- 
sults obtained with indicated parametrizations of the RMF 
Lagrangian in the Ar isotopes. The structure of the ground 
states is shown by the Nilsson labels. Only in the case when 
the structure of the ground state is the same as in the NL3 
parametrization, the results obtained with other parametriza- 
tions are shown. 

of different (quasi)particle states in model calculations is 
also needed. 

While there was a considerable interest in the study of 
time-odd mean fields in odd- and odd-odd mass nuclei at 
no rotation within the Skyrme EDF [24|, |25[, relatively 
little is known about their role in the framework of the 
CDFT. So far the impact of time-odd mean fields on 
binding energies has been studied in the CDFT frame- 
work only in odd-mass nuclei around doubly magic spher- 
ical nuclei in Ref. 15 ill, and in few deformed nuclei around 
32 S [III and 254 No|26j]. Only last two publications treat 
time-odd mean fields in a fully self-consistent manner; 
the breaking of Kramer's degeneracy is neglected in Ref. 



A. Binding energies in light nuclei 

The impact of NM on the binding energies of light 
odd-mass nuclei is shown in Fig. [T] One can see that in 
all cases the presence of NM leads to additional binding, 
the magnitude of which is nucleus and state dependent. 
The absolute value of this additional binding is typically 
below 200 keV and only in some lower mass nuclei it 
reaches 300 keV. On the average, the magnitude of ad- 
ditional binding due to NM correlates with the mass of 
nucleus; it is largest in lightest nuclei and smallest in the 
heaviest nuclei. For each isotope chain, it is largest in the 
vicinity of the proton-drip line and smallest in the vicin- 
ity of neutron-drip line. Fig. [2] shows that it only weakly 
depends on the RMF parametrization; this is also seen 
in the analysis of terminating states in Ref. [l^| . In both 
cases, the biggest deviation from the NL3 results is ob- 



served in the case of the NLSH parametrization. 

It is interesting to compare these results with the ones 
obtained in the Skyrme EDF (see Fig. 4 in Ref. [IJ]). 
The modifications of total binding energy due to time- 
odd mean fields are given by the E quantity in Ref. 
[U, which is an analog of the E NM - E WNM quantity. 
The general dependence of both quantities on N — Z 
is similar in odd-mass nuclei apart of few cases such as 
43 Ti and 43 Sc in SLy4 Skyrme EDF (Fig. 4 in Ref. [H). 
Neither of them indicates the enhancement of time-odd 
mean fields in the vicinity of the N = Z line. In addition, 
the absolute values of E to and E NM - E WNM quantities 
are similar being below 300 keV in the majority of the 
cases. The principal difference between the RMF and the 
Skyrme EDF lies in the fact that time-odd mean fields 
are always attractive and show very small dependence on 
the parametrization in the RMF calculations (this is also 
supported by the analysis of terminating states, see Ref. 
[19(), while they can be both attractive (SLy4 force) or 
repulsive (SIII force) and show considerable dependence 
on the parametrization in Skyrme EDF (Ref. [24]]). 



B. Binding energies in the Ce (Z = 58) isotopes 

The role of time-odd mean fields is studied here in 
medium mass Ce isotopes in order to facilitate the com- 
parison with the results obtained within the Skyrme EDF 
with the SLy4 force in Ref. [25[ . This reference represents 
the most detailed study of time-odd mean fields in odd- 
mass nuclei within the Skyrme EDF. We consider lowest 
one-particle states of positive and negative parities, while 
Ref. [2 51 ] studies only lowest one-particle state in each nu- 
cleus. 

Figs. [3] and [4] show the additional binding due to NM. 
The comparison with the Skyrme EDF results of Ref. 
[IH reveals the number of important differences. First, 
similar to the results in light nuclei (Sect. IIII Ap and in 
actinide region (Sect. VI H in Ref. [26|), time-odd mean 
fields are attractive in the RMF calculations for the Ce 
isotopes. On the contrary, they are repulsive in the SLy4 
parametrization of the Skyrme EDF [25J . Note that the 
SLy4 force produces attractive time-odd mean fields in 
light nuclei (Ref. [24[). This mass dependence of the ef- 
fects of time-odd mean fields in the Skyrme EDF may be 
due to the competition between isovector and isoscalar ef- 
fects [25{ . The average absolute magnitude of the change 
of binding due to time-odd mean fields in the RMF calcu- 
lations is only half of the one seen in the Skyrme calcula- 
tions with the SLy4 parametrization. It was also checked 
on some examples that additional binding due to NM 
only weakly depends on the parametrization of the RMF 
Lagrangian. 

Second, the results of calculations do not reveal strong 
dependence of additional binding due to NM on defor- 
mation. For example, the deformation of the ^[615]ff/2 
state in the 173_181 Ce chain changes drastically from 
P 2 ~ 0.23 down to /3 2 ~ 0.06 (Fig. OH bottom panel), but 
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FIG. 3: The impact of NM on binding energies of positive 
parity states in odd mass Ce (Z = 58) nuclei. Upper panel 
shows the additional binding E NM ~ E WNM due to NM and 
its configuration dependence. The structure of the states is 
shown by the Nilsson labels; the states at and to the right 
of the Nilsson label up to the next Nilsson label have the 
same structure. The bottom panel shows the corresponding 
deformations of the single-particle configurations. The calcu- 
lations have been performed with the NL3 parametrization of 
the RMF Lagrangian. They cover the nuclei from proton-drip 
line up to neutron-drip line. 



the additional binding due to NM remains almost the 
same (Fig. H top panel). The i/[523]7/2 and z/[505]ll/2 
states are another examples of this feature (Fig. [J). 

Third, the binding energy modifications due time-odd 
mean fields are completely different in the RMF and 
Skyrme EDF calculations. In the Skyrme EDF calcu- 
lations, the magnitude of these modifications correlates 
with three properties of the blocked orbital, namely, in 
decreasing order of importance they are: a small Q quan- 
tum number, a down-sloping behavior of the energy of the 
single-particle state with mass number A, and a large 
total angular momentum j for the spherical shell from 
which the single-particle state originates [25| . On the 
contrary, the RMF calculations do not reveal this type 
of correlations. Indeed, the single-particle states which 
have the largest changes in binding energies due to NM 
(\E NM - E WNM \ > 0.1 MeV) are [413]5/2, [404] 7/2, 
[640] 1/2, [631]3/2, [505] 9/2 and [501] 1/2. 

Neutron current distributions in one-neutron configu- 




120 130 140 150 160 170 180 

Mass number A 

FIG. 4: The same as in Fig. [3] but for negative parity states. 



rations of selected nuclei having similar quadrupole de- 
formations are shown in Figs. O and [5] They are predom- 
inantly defined by the currents generated by the blocked 
orbitals and are characterized by the complicated pat- 
terns in different cross-sections of nucleus. However, 
there are clear correlations between the patterns of the 
currents in the y — z plane (z is the symmetry axis and 
x is the rotation axis in the CRMF) and the projection 
VLu of the angular momentum of the blocked orbital on 
the axis of symmetry (Fig. [5]). At fi^ = 1/2, the currents 
show circulations (vortices) which are concentrated in the 
central region of nucleus. However, with increasing Q&j, 
the currents are pushed away from the axis of symme- 
try of the nucleus towards the surface area. In addition, 
the strength of currents correlates with flu- as follows 
from the values of factor F the strongest currents appear 
for low fif,;, while the configurations with large Flu are 
characterized by weak currents. In the x — y plane, the 
majority of the configurations show the current pattern 
(although with different strengths of the currents) visible 
on Fig. [BJd, while in the x — z plane the currents are typ- 
ically directed towards (away) from the z = axis (see, 
for example, Fig. H^i). Figs.[S]i,e show that the change of 
the signature of the blocked orbital leads to a change in 
the direction of the currents. 
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FIG. 5: (Color online) Neutron current distributions j n {r) in intrinsic frame in the y — z plane at x — 0.48 fm for several 
one- neutron configurations in the Ce nuclei. The Nilsson quantum numbers [Nn z A]Qbi indicate the blocked state; the sign 
of its signature r is shown by the superscript. The currents in panels (d) and (e) are plotted at arbitrary units for better 
visualization. The currents in other panels are normalized to the currents in panels (d) and (e) by using factor F. This factor 
is chosen in such way that the current distribution for every nucleus is clearly seen. The shape and size of the nucleus are 
indicated by density lines which are plotted in the range 0.01 - 0.06 far 3 in step of 0.01 fm" 3 . The panels are arranged in 
such way that the Qm value of the Nilsson label of the blocked state increases on going from panel (a) to panel (h) . 



C. Mass and particle number dependences of 
additional binding due to NM 

Mass, neutron and proton number dependences of ad- 
ditional binding due to NM (the \E NM - E WNM \ quan- 
tity) are presented in Figs. [7] and [8] These figures 
are based on the results obtained in Sects. IIII Al and 
IIIIBI and on some extra calculations for odd-mass nu- 
clei around 249 Cf. The latter calculations were per- 
formed in order to check the impact of pairing on the 
\E NM ~ E WNM \ quantity and to add some high- A, Z, N 
data to better define particle number dependences of 
this quantity. These calculations without pairing were 



compared with the Relativistic Hartree-Bogoliubov re- 
sults for some one-quasiparticle states shown in Table 
IV of Ref. [26}. Although the pairing decreases addi- 
tional binding due to NM in most of the cases, there 
are still one- (quasi) particle configurations in which the 
\E NM - E WNM I quantity is smaller in the calculations 
without pairing. This is a consequence of the compli- 
cated nature of the E NM - E WNM quantity defined by 
(i) the interplay of time-odd mean fields and the polar- 
ization effects (Sect. HV]) and by (ii) the differences in the 
impact of pairing on different terms of total energy. 

The calculated \E NM - E WNM \ quantities were fitted 
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by simple parametrization 



AE = 




where Q is equal either to proton Z , neutron iV or mass 
A numbers. Note that the \E NM - E WNM \ values from 
odd-proton (odd-neutron) nuclei were used in the fit of 
Z — (N-) dependence of AE, while all available values 
were used in the fit of A— dependence. The results of the 
fits are shown by solid lines in Figs. and [S] One can 
see that the powers a are similar for different fits (proton, 
neutron or mass) and it is likely that some differences be- 
tween them are due to limited set of data used in the fit. 
On the other hand, the magnitudes c differ considerably 
between proton and neutron quantities. More extensive 
calculated data for \E NM ~ e wnm \ will definitely im- 
prove the accuracy of such fits and may allow to define 
the sources of observed difference in the magnitudes c 
between proton and neutron quantities. However, such 
study is beyond the scope of the current manuscript. 
The three-point indicator [54[ 

A^(N) = 2£ [B(N - 1) + B(N + 1) - 2B(N)} (26) 

is frequently used to quantify the odd-even staggering 
(OES) of binding energies. Here tin = (— 1) is the num- 
ber parity and B(N) is the (negative) binding energy of a 
system with N particles. In Eq. (|26p. the number of pro- 
tons Z is fixed, and N denotes the number of neutrons, 
i.e. this indicator gives the neutron OES. The factor de- 
pending on the number parity 7Tjv is chosen so that the 
OES centered on even and odd neutron number N will 
both be positive. An analogous proton OES indicator 
A^\Z) is obtained by fixing the neutron number N and 
replacing N by Z in Eq. fl2£§. 

The A^(N) (and similarly A<- 3 \Z)) quantity will be 
modified in the presence of time-odd mean fields as 

A^ (N) = A$ TO (N) + SE TO (27) 

where the subscripts 'TO' and 'WTO' indicate the values 
obtained in the calculations with and without time-odd 
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FIG. 7: (Color online) Mass dependence of additional binding 
due to NM. Open and solid circles are used for odd-proton and 
odd-neutron nuclei, respectively. 



mean fields, and SEto is the contribution coming from 
time-odd mean fields. If the A^ 3 '(7V) quantity is cen- 
tered at odd-iV nucleus, the SEto quantity represents 
the change of binging energy of this odd-mass nucleus 
induced by time-odd mean fields. This is because the 
time-odd mean fields have no effect on binding energies 
of the ground states of even-even nuclei. Note that with 
such selection SEto is negative if time-odd mean fields 
provide additional binding in odd-mass nucleus. 

In the CDF theory, the SEto quantity is equal to 
E — E WNM and thus it is always negative: this result 
does not depend on the RMF parametrization (see Sect. 
UTO] for the dependence of the E NM - E WNM quantity 
on the RMF parametrization). In addition, the magni- 
tude of the SEto quantity depends only weakly on the 
RMF parametrization. On the contrary, the sign and the 
magnitude of SEto depends strongly on the parametriza- 
tion in the Skyrme EDF calculations. For example, in 
the calculations with the SLy4 force the SEto quantity 
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is positive for medium mass nuclei (Refs. [Ill [H, Ull ) but 
negative in light nuclei (Ref. O n the other hand, 

the SEto quantity will be positive in light nuclei in the 
calculations with the SIII parametrizations [24| . 

It is interesting to compare the averaged effects of time- 
odd mean fields as given by the AE quantity with the 
experimental global trends for OES as shown by dashed 
lines in Fig. 2 in Ref. [55J]. The latter trends were ob- 
tained using phenomenological parametrization with the 
same functional dependence as in Eq. ([25| with c = 4.66 
MeV (4.31 MeV) and a = 0.31 for neutron (proton) data 
sets. The comparison of theory and experiment suggests 
that time-odd mean field contributions into OES can be 
as large as 10% in light systems and around 5-6% in 
heavy systems. These are non-negligible contributions 
which have to be taken into account when the strength 
of pairing interaction is defined from fits to experimental 
OES. The analysis of the Sn isotopes in Ref. [29| showed 
that the time-even and time-odd polarization effects in- 
duced by the odd nucleon produce about 30% reduced 
OES as compared to the standard spherical calculations. 
As a consequence, an enhancement of pairing strength 
by about 20% is required to compensate for that effect. 
Our calculations show much smaller reduction of OES in 
part because the polarization effects in time-even channel 
are already taken into account in the calculations with- 
out NM. In addition, the polarization effects in time-odd 
channel (in particular, the brea king of Kramer's degen- 
eracy which is neglected in Ref. 29]) are more correctly 
treated in the current manuscript (see Sect. HV|) . Thus, 
the current calculations suggest that much smaller in- 
crease of the strength of pairing (by approximately 5%) 
would be required to compensate the reduction of OES 
due to time-odd mean fields. 



Detailed analysis of the single-particle spectra in 119 Ce 
and 123 Xe reveals general features which are also found 
in other nuclei. The 119 Ce nucleus is axially symmetric 
(7 = 0°) while 123 Xe is triaxial with 7 = -26°. This 
difference in the symmetry of nucleus results in important 
consequences: the energy splittings appear in all single- 
particle states in triaxial nuclei, while only the states 
with £1 = flu (fi is the projection of the total angular 
momentum of the particle on the axis of symmetry and 
the subscript 'bl ' indicates the blocked state) experience 
such splittings in axially symmetric nuclei. The former 
feature is due to the fact that Q is not good quantum 
number in triaxial nuclei and each single-particle state 
represents a mixture of the basic states with different 
values of f2. 

It is important to mention that the occupied and unoc- 
cupied states as well as proton and neutron states show 
energy splittings (Fig. [9]). The splittings of the proton 
and neutron states of the same structure are similar. This 
is because the largest contribution to magnetic potential 
(Eq. ([5])) is due to space- like components of the w-meson 
fields which do not depend on the isospin. In addition, 
the occupied state is always more bound than its unoc- 
cupied time-reversal counterpart. 

The change of the signature of the blocked state leads 
to the inversion of the signatures in other pairs of time- 
reversal orbitals (compare columns (1) and (3) in Fig. [5]). 
This is due to the change of the direction of the currents 
caused by the change of the signature of the blocked state 
(compare Fig. [SJi with Fig. [5J3). However, the additional 
binding due to NM (the E NM - E WNM quantity) does 
not depend on the signature of the blocked state in odd- 
mass nuclei. 



IV. THE MECHANISM OF ADDITIONAL 
BINDING DUE TO NM IN ODD-MASS NUCLEI 

In the current section a detailed analysis of the impact 
of NM on the energies of the single-particle states and on 
different terms in the total energy expression (Eq. (fT5|l ) is 
performed in order to better understand the microscopic 
mechanism of additional binding due to NM. We use the 
z/5/2[413] configuration of 119 Ce as an example in this 
analysis. 



A. Energy splittings of time-reversal counterpart 
single-particle states in the presence of NM 

Fig. [5] shows that the presence of time-odd mean fields 
leads to the energy splitting AE sp i it (i) of the single- 
particle states which are time-reversal counterparts. This 
corresponds to the removal of the Kramer's degener- 
acy of these states. One of these states moves up by 
~ AE sp ut/2 as compared with its position in the absence 
of NM, while another down by w AE sp ut/2. 



B. Polarization effects induced by NM 

The polarization effects induced by NM are investi- 
gated by considering its impact on different terms of the 
total energy (Eq. (fT5|) ). The results of this study are 
shown in Table [TTJ One can see that the total energy 
terms can be split into two groups dependent on how 
they are affected by NM. The first group includes the 
Eunl, Ej, Ep and EqouI terms which are only weekly 
influenced by NM, and thus, they will not be discussed 
in detail. 

The second group is represented by the E part , E a , E^ 
and E^ terms which are strongly affected by NM. The 
£"5 term is directly connected with the nucleonic currents 
(see Eq. (|2"TjO . The E a and terms depend only in- 
directly on the time-odd mean fields: the minimization 
of the total energy in the presence of time-odd terms 
leads to a very small change of equilibrium deformation 
induced by NM. The quadrupole and hexadecapole mo- 
ments change by 10~ 4 of their absolute value when the 
NM is switched on; a similar magnitude of changes is seen 
also in E a and E^ . One should keep in mind that only 
the E a + E^ quantity has a deep physical meaning since 
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FIG. 8: (Color online) The same as in Fig. [7] but for neutron and proton dependences of additional binding due to NM. 
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FIG. 9: (Columns (1) and (3)) The energy splittings AE sp u t 
between different signatures of the single-particle states in the 
presence of NM. The results of the calculations are shown for 
the configurations of 119 Ce in which either the ^[413]5/2~ 
(column (1)) or ^[413]5/2 + (column (3)) states are blocked. 
These signatures are degenerated in energy in the calculations 
without NM (column (2)). Note that the single-particle states 
of interest are shown at arbitrary absolute energy in the col- 
umn (2). Solid (open) circles indicate occupied (unoccupied) 
states. Solid (dotted) lines are used for the r — +i (r = — i) 
states. 



it defines a nucleonic potential; this sum is modified by 
NM only on -177 keV. 

The largest modification (by —410 keV) is seen in the 
Epart energy, with the half of it coming from the change 
of the single-particle energy (by « —200 keV) of the 
blocked orbital (the i/[413]5/2 orbital) in the presence 
of NM. Note that since both signatures of other pairs 
of time-reversal orbitals below the Fermi level are occu- 
pied, the large energy splittings E sp n t seen for some of 



TABLE II: The impact of NM on different terms of the to- 
tal energy (Eq. ([15)1 ) in the [413]5/2 configuration of 119 Ce. 
Second column shows the absolute energies [in MeV] of dif- 
ferent energy terms in the case when NM is neglected. Third 
column shows the changes [in MeV] in the energies of these 
terms induced by NM. 



Quantity 


E t 


tt^NM r?WNM 
Ei — ti>i 


1 


2 


3 


Epart 


-2849.889 


-0.410 




-17079.532 


-2.231 


EcrNL 


343.341 


-0.017 


El 


14356.156 


2.054 


El 


0.0 


-0.124 


Ej 


2.044 


0.003 


E S P 


0.0 


-0.010 


Ecoui 


481.196 


0.017 


E cm 


-6.252 


0.0 


Etot 


-959.349 


-0.104 



them do not have a considerable impact on E part (see 
Eq. (|16p ) since this splitting is nearly symmetric with 
respect to the position of these orbitals in the absence 
of NM. Thus, the rest of the modification of E part is re- 
lated to small changes in the single-particle energies of 
occupied states caused by the changes in the equilibrium 
deformation induced by NM. 

This detailed analysis clearly indicates that the E NM — 
E WNM quantity is defined by both time-odd mean fields 
and the polarization effects in time-even mean fields in- 
duced by time-odd mean fields. E NM - E WNM = -104 
keV is a result of near cancellation of the contributions 
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FIG. 10: The ratio AE spUt /(E NM - E WNM ) in the Ce iso- " 6 °0 5 10 15 

topes. The structure of the blocked states is shown by the Radius r [fm] 

Nilsson labels; the states at and to the right of the Nilsson 
label up to the next Nilsson label have the same structure. 



due to fermionic (—410 keV) and bosonic (—306 keV) 
degrees of freedom. Note that the latter appears with 
negative sign in Eq. (|15[) . The fermionic degrees of 
freedom are represented by the E part and E cm terms, 
while other terms of the total energy are related to the 
bosonic degrees of freedom. The fermionic contribution 
into E NM — E WNM is defined by more or less equal con- 
tributions from time-odd mean fields and polarization 
effects in time-even fields. On the contrary, the time-odd 
mean fields define only » 1/3 {E% = -0.124 keV) of the 
bosonic contribution into E NM - E WNM , while the rest 
is due to polarization effects in time-even mean fields. 

It turns out that these contributions are highly cor- 
related as can be seen from the ratio AE sp u t /(E NM — 
E WNM ) in the Ce isotope chain (Fig. US). AE spUt de- 
pends only on time-odd mean fields in fermionic chan- 
nel, while (E NM - E WNM ) depends both on time-odd 
mean fields and polarizations effects in time-even mean 
fields in fermionic and bosonic channels. One can see 
that AE spHt /(E NM -E WNM ) w 4 for majority of nuclei. 
Similar relation exists also in the Skyrme EDF calcula- 
tions for the Ce isotopes (see Eq. (7) in Ref. (25j|). 



V. THE IMPACT OF TIME-ODD MEAN 
FIELDS ON THE PROPERTIES OF 
PROTON-UNSTABLE NUCLEI 

The blocked state has always lower energy than its 
unoccupied time-reversal counterpart in the calculations 
with NM; this fact does not depend on the signature of 
the blocked state (Sect. \TV Aj) . The energy of blocked 
state in the presence of NM is lower by « AE sp u t /2 than 
the energy of the same state in the absence of NM. This 
additional binding will affect the properties of the nuclei 
in the vicinity of the proton-drip line via two mechanisms 
discussed below. They are schematically illustrated in 



FIG. 11: Schematic illustration of the impact of time-odd 
mean fields on the properties of odd-proton nuclei in the vicin- 
ity of proton drip line. The single-proton states, involved 
in the mechanisms discussed in the text, and proton nucle- 
onic potential (which also includes the Coulomb potential) 
are shown in the figure. 



Fig. HQ 

In the first mechanism, the nucleus, which is proton 
unbound (state A in Fig. [H]) in the calculations with- 
out NM, becomes proton bound in the calculations with 
NM (state A' in Fig. [TTj) . The necessary condition for 
this mechanism to be active is the requirement that the 
energy of the single-particle state in the absence of NM 
is less than AE sp n t /2. This mechanism can be active 
both in the ground and excited states of the nuclei in the 
vicinity of the proton-drip line. 

In the second mechanism, the energy of the single- 
particle state (state B' in Fig. [IT]) is lower in the presence 
of NM, but the state still remains unbound. This will 
affect the decay properties of proton emitters and the 
possibilities of their observation. Indeed, the lowering 
of the energy of the single-proton state will decrease the 
probability of emission of the proton through combined 
Coulomb and centrifugal barrier. Many results of the 
physics of proton emitters are conventionally expressed 
in terms of the Q p energies which depend on the dif- 
ference of binding energies of parent (odd-proton) and 
daughter (even-proton) nuclei. Note that for simplicity 
we consider here only even-TV nuclei. NM leads to an 
additional binding in odd-proton nucleus but it does not 
affect the binding of even-proton nucleus. Thus, the Q p 
values are lower by the value of this additional binding 
when the NM is taken into account. 

Two consequences follow from lower Q p values. First, 
experimental observation of proton emission from the nu- 
cleus will become impossible if the Q p value moves out- 



12 



side the Q p window favorable for the observation of pro- 
ton emission or becomes possible if the Q p value moves 
into the Q p window favorable for the observation of pro- 
ton emission. The size of the Q p window for rare-earth 
proton emitters is about 0.8 — 1.7 MeV, while it is much 
smaller in lighter nuclei [EH . Large Q p values outside 
this window result in extremely short proton-emission 
half- lives, which are difficult to observe experimentally. 
On the other hand, the decay width is dominated by 
f3 + decay for low Q p values below the Q p window. This 
consequence of the lowering of Q p due to NM is espe- 
cially important in light nuclei where the impact of NM 
on binding energies is especially pronounced and the Q p 
window is narrow. 

Second, the lowering of the Q p values due to NM will 
increase the half-lives of proton emitters. For example, 
the lowering of Q p due to NM will be around 50 keV 
in rare-earth region since this is typical value of addi- 
tional binding due to NM in odd-mass nuclei of this re- 
gion (Sect. IhTB]) . This can increase the half-lives of pro- 
ton emitters by a factor of w 2 at the upper end of the 
Q p window and by a factor of w 4 at the bottom end 
of the Q p window (see Fig. 5 in Ref. [58(). The effects 
of NM have been neglected in the existing RHB studies 
of p roton emitters with Z > 50 (see, for example, Ref. 
[59() but this should not introduce significant error in 
this mass region. 

On the other hand, the impact of NM can be dramatic 
on the half-lives of proton emitters in lighter nuclei. This 
is due to two factors, namely, (i) the general increase 
of additional binding due to NM and the magnitude of 
AE sp iit with decreasing mass and (ii) the narrowing of 
the Q p window with the decrease of mass due to the 
lowering of the Coulomb barrier. This can be illustrated 
by several examples. The change in proton energy of 
around 300 keV in 69 Br causes a change in the proton 
decay lifetime of 11 orders of magnitude - This effect 
is even more pronounced in lighter systems. The half-life 
window of 10 to 10~ 4 s corresponds to proton energies 
of 100 - 150 keV in nuclei around Z = 20 [60], while 
the variation of the Q p value between 3 to 50 keV in 
7 B changes the half-lifes by 30 orders of magnitude [58| . 
The energy changes quoted in these examples are either 
of similar magnitude or even smaller as compared with 
the changes of the energies of single-proton states and the 
Q p values induced by NM. As a result, one can conclude 
that the effects of time-odd mean fields have to be taken 
into account when attempting to describe the properties 
of proton emitters in light nuclei. 



VI. ODD-ODD MASS NUCLEI: A MODEL 
STUDY OF IMPACT OF NUCLEAR 
MAGNETISM ON BINDING ENERGIES. 

The nuclei around 32 S in superdeformed minimum are 
considered in the present Section. Their selection is 
guided in part by the desire to compare the CRMF re- 
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FIG. 12: Neutron single-particle energies (routhians) as a 
function of the rotational frequency Q x . They are given 
along the deformation path of the doubly-magic SD config- 
uration 7r3 2 i^3 2 in 32 S. Solid, short-dashed, dot-dashed and 
dotted lines indicate (n — +, r = — i), (n = +, r — 
(-7T = — , r = +i) and [it = — , r = —i) orbitals, respectively. 
At Q x = 0.0 MeV, the single-particle orbitals are labeled by 
means of the asymptotic quantum numbers [Nn z A]fi (Nilsson 
quantum numbers) of the dominant component of the wave 
function. 



suits with the ones obtained in the Skyrmc EDF in Ref. 
[6l[ , where the signature separation induced by time-odd 
mean fields has been found in the excited SD bands of 
32 S. The CRMF calculations have been performed for 
some SD configurations in 32 S and in neighboring nuclei. 
The starting point is the doubly magic SD configura- 
tion tt3 2 z/3 2 in 32 S (further 'SD core') (see Ref. (HI) in 
which all single-particle orbitals below the N = Z = 16 
SD shell gaps are occupied (Fig. [T2"|) . Then the config- 
urations in the nuclei under consideration (Fig. I13p are 
created by either adding particles into the [202]5/2 ± or- 
bitals) or/and creating holes in the [330]l/2 ± orbitals: 
these are the orbitals active in signature-separated con- 
figurations discussed in Ref. (6l| . 

Similar to the results shown in Sects. IIII Al and IIII B[ 
the NM leads to additional binding in one-particle con- 
figurations of odd mass nuclei (the configurations in 3 S 
and 33 CI created by adding a particle to the SD core or 
the configurations in 31 P and 31 S created by removing a 
particle from the SD core, see Fig. [T3"|) . This additional 
binding does not depend on the signature of the single- 
particle state. 

Fig. EH shows that additional binding due to NM is 
smaller for the proton states as compared with the neu- 
tron ones. For example, the single-particle configurations 
in 31 P and 31 S are built on the same Nilsson state. How- 
ever, additional binding due to NM is smaller for the pro- 
ton state ( 31 P) than for the neutron state ( 31 S). Similar 
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FIG. 13: The impact of NM on binding energies of different 
single-particle configurations under study. The energies of 
single-particle configurations calculated without NM are nor- 
malized to zero. Short lines show the magnitude of additional 
gain (negative energies) or loss (positive energies) in binding 
energies in the presence of NM. The configurations are labeled 
by the particle (p) and/or hole (h) states with respect to the 
32 S SD core and the total signature r of the configuration. 



FIG. 14: (Color online) The impact of NM on binding ener- 
gies of the lowest two-particle configurations in odd-odd Al 
nuclei. Upper panel shows the E NM — E WNM quantity for 
different signatures. The structure of the states are shown by 
the Nilsson labels only in the cases when two-particle config- 
urations are near-prolate. The same state is blocked in the 
proton subsystem of all nuclei. The bottom panel shows the 
fh and 7-deformations of nuclei in the configurations under 
study. 



situation exists also in 33 S and 33 CI. These results are 
consistent with a general systematics (Sect. MI CP which 
shows that additional binding due to NM is smaller in 
the proton subsystem than in the neutron one. The re- 
distribution of the proton density due to Coulomb inter- 
action may be one of possible reasons for such difference 
between proton and neutron states. 

The situation is more complicated in odd-odd nuclei 
( 30 P and 34 C1) in which considerable energy splitting be- 
tween the r = +1 and r = — 1 configurations is obtained 
in the calculations. The microscopic mechanism of bind- 
ing modifications is illustrated in Table [TTTI on the exam- 
ple of configurations A and B in 34 C1. 

NM provides additional binding of around 0.4 MeV in 
the configuration A which has signature r = — 1. The 
proton and neutron currents due to the occupation of 
proton and neutron 5/2 [202] ~ states are in the same 
direction which results in appreciable baryonic current. 
This baryonic current leads to sizable modifications in 
the E part , E a , E%, El terms (Table Hill). These are pre- 
cisely the same terms which are strongly affected by NM 
in odd-mass nuclei, see Sect. II VB] The fermionic contri- 
bution into E NM - E WNM (the E part term) is denned 
by more or less equal contributions from time-odd mean 
fields and polarization effects in time-even mean fields. 
On the contrary, the time-odd mean fields define only 
w 1/3 (E u = -0.413 MeV) of the bosonic contribution 
into E NM — E WNM , while the rest is due to polarization 
effects in time-even mean fields (the E a , E£ terms). 

On the contrary, the NM leads to the loss of binding in 
the configuration B which has r = +1. In this configura- 
tion, the proton and neutron currents due to 7r[202]5/2 + 
and j/[202]5/2~ states are in opposite directions, so the 




N-Z 



FIG. 15: (Color online) The same as in Fig. [14] but for the 
lowest two-particle configurations in odd-odd CI nuclei. 



total baryonic current is very close to zero. As a result, 
the impact of NM is close to zero for the majority of the 
terms in Eq. (fTS]) (see Table [TTT]) . The only exception 
is the Ep term which represents space-like component of 
the isovector-vector p-field. This term depends on the dif- 
ference of proton and neutron currents (Eq. (|22p ). which 
for the present case of opposite currents gives non-zero 
result. As follows from Table ITTT] this term is predomi- 
nantly responsible for the loss of binding due to NM in 
the configuration B. 

Fig. Q2] also shows the results for the 4-particle excited 
SD states Tr(a+b)+i , (a+b) in 32 S, for which the calculated 
rotational structures display the signature separation in- 
duced by time-odd mean fields [52|, [6l| . The configura- 
tions are formed by exciting proton and neutron from 
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TABLE III: The changes in different terms of the total energy (Eq. (jT5j ) in the SD configurations A = 7r[202]5/2 <g> 
//[202]5/2"(r = -1) and B = tt[202]5/2 + <g> i/[202]5/2 _ (r = +1) of 34 C1 induced by NM (third and fourth columns). Second 
column shows the absolute energies [in MeV] of different energy terms in the case when NM is neglected. The configurations 
are given with respect of doubly magic SD configuration in 32 S. 



Quantity 


E^ NM (A,B) 


Ei fM (A)-Ej VNM (A) 


Ef fM (B)-Ef NM (B) 


1 


2 


3 


4 


Epart 


-836.044 


-1.47 


-0.004 


E a 


-4416.351 


-7.824 


+0.004 


EoNL 


84.864 


-0.036 


-0.001 


El 


3698.757 


7.126 


-0.003 


El 


0.0 


-0.413 


0.0 


E T P 


0.061 


0.0 


0.0 


E S P 


0.0 


0.0 


-0.043 


Ecoui 


59.839 


0.051 


-0.002 


R 


-9.492 


0.0 


0.0 


Etot 


-272.705 


-0.374 


0.041 
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FIG. 16: (Color online) The impact of NM on binding en- 
ergies of two-particle states in odd-odd nuclei in the vicinity 
of the N = Z line as a function of proton number Z. Panel 
(a) shows the results for two-particle configurations with dif- 
ferent proton and neutron states, while panel (b) shows the 
results for two-particle configurations with the same proton 
and neutron states. The structure of the states are shown by 
the Nilsson labels only in the cases when two-particle config- 
urations are near-prolate. Note that in panel (b) only one 
Nilsson label is shown since proton and neutron states have 
the same structure. The results are shown only in the cases 
when the convergence has been achieved for both N = Z and 
N = Z - 2 (or N = Z + 2) nuclei. 



the [330] 1/2" orbitals below the N = 16 and Z = 16 
SD shell gaps into the [202] 5/2 ± orbitals located above 
these gaps. They have the 7r3 1 ^3 1 structure in terms of 
intruder orbitals. When NM is neglected these four con- 
figurations are degenerated in energy. This degeneracy 



is broken and additional binding, which depends on the 
total signature of the configuration (0.907 MeV for the 
r = +1 configurations and 0.468 MeV for the r = — 1 con- 
figurations in the calculations with the NL3 parametriza- 
tion), is obtained when NM is taken into account. The 
NL1 and NLSH parametrizations of the RMF Lagrangian 
give very similar values of additional binding due to NM. 
The essential difference between the relativistic and non- 
relativistic calculations lies in (i) the size of the energy 
gap between the r = +1 and r — — 1 configurations and 
(ii) the impact of time-odd mean fields on the energy of 
the r = — 1 states. This energy gap is about 2 MeV in the 
Skyrme EDF calculations with the SLy4 force [6l| , while 
it is much smaller being around 0.45 MeV in the CRMF 
calculations with the NL1, NL3 and NLSH parametriza- 
tions. The energies of the r = — 1 states are not affected 
by time-odd mean fields in the Skyrme EDF calculations 
[6l| , while appreciable additional binding is generated by 
NM for these states in the CRMF calculations (Fig. [H]). 

Figs . [T4l and [T5l show the results of the calculations for 
the ground state two-particle configurations in odd-odd 
Al and CI nuclei. The calculations suggest that signa- 
ture separation due to time-odd mean fields is expected 
also in two-particle configurations at zero or low exci- 
tation energies in odd-odd nuclei. The signature sep- 
aration is especially pronounced in the N — Z 26 Al 
(the 7r[202]5/2®i/[202]5/2 configuration) and 34 C1 nuclei. 
This is because the proton and neutron currents in these 
configurations are almost the same both in strength and 
in spatial distribution. As a result, their contribution to 
total energy is large in the case when these currents are 
in the same direction (the r = — 1 configurations) and 
close to zero when these currents are in opposite direc- 
tions (the r = +1 configurations). Note that 26 Al is axi- 
ally deformed while the CI is triaxially deformed with 



15 



7 ~ 30°. However, both of them show the enhancement 
of signature separation at N = Z. 

The signature separation is rather small for the ma- 
jority of nuclei away from the N = Z line. This is a 
consequence of the fact that the strength of the currents 
in one subsystem (and thus the impact of NM on binding 
energies) is much stronger that in another subsystem. As 
a result, there is no big difference (large signature sepa- 
ration) between the cases in which the proton and neu- 
tron currents are in the same and opposite directions. 
However, some nuclei away from the N = Z line also 
show appreciable signature separation. These are 38 Al 
and 38 . 48 . 50 CI nuclei (Figs. ITS IT5]) for which the strengths 
of proton and neutron currents (but not necessarily the 
spatial distribution of the currents) are of the same order 
of magnitude. 

It was suggested in Ref. [24| that the effects of time-odd 
mean fields are enhanced at the N = Z line. However, 
Fig. [TB] clearly shows that the enhancement of signature 
separation is not restricted to the N = Z line. Indeed, 
signature separation of the same two-particle configura- 
tions are very similar in the N = Z and N = Z ± 2 nuclei 
despite the fact that sometimes the deformations of com- 
pared nuclei differ appreciably. There are considerable 
signature separation in two-particle configurations based 
on the same proton and neutron states in the N = Z 
and N = Z — 2 nuclei (Fig.fTBb). On the other hand, al- 
most no signature splitting is observed in the N = Z and 
N = Z + 2 nuclei when two-particle configurations are 
based on different proton and neutron states (Fig. ITBa). 
This suggests that the enhancement of signature splitting- 
is due to similar proton and neutron current distributions 
(see discussion in previous paragraph) . 

When considering odd-odd nuclei one has to keep in 
mind that the present approach takes into account only 
the part of correlations between blocked proton and neu- 
tron and neglects the pairing. In particular, the residual 
interaction of unpaired proton and neutron leading to 
the Gallagher-Moshkowski doublets of two-quasi par ticle 
states with K > = Q p + Q n and K < — \fl p — Q n | |62l . [63| 
is not taken into account. Thus, future development of 
the model is required in order to compare directly the 
experimental data on odd-odd nuclei with calculations. 
This question will be discussed in more details in a forth- 
coming manuscript [3fJ. 



VII. CONCLUSIONS 

Time-odd mean fields (nuclear magnetism) have been 
studied at no rotation in a systematic way within the 
framework of covariant density functional theory. The 
main results can be summarized as follows: 

• In odd-mass nuclei, nuclear magnetism always leads 
to an additional binding indicating its attractive 
nature in the CDFT. This additional binding only 
weakly depends on the parametrization of the RMF 



Lagrangian. On the contrary, time-odd mean fields 
in Skyrme EDF can be attractive and repulsive and 
show considerable dependence on the parametriza- 
tion of density functional. This additional binding 
is larger in odd-neutron states than in odd-proton 
ones in the CDFT framework. The underlying mi- 
croscopic mechanism of additional binding due to 
NM has been studied in detail. 

• Additional binding due to NM can have a profound 
effect on the properties of odd-proton nuclei in the 
ground and excited states in the vicinity of the 
proton-drip line. In some cases it can transform the 
nucleus which is proton unbound (in the calcula- 
tions without NM) into the nucleus which is proton 
bound. This additional binding can significantly af- 
fect the decay properties of proton unbound nuclei 
by (i) increasing the half-lives of proton emitters 
(by many orders of magnitude in light nuclei) or 
(ii) moving the Q p value inside or outside the Q p 
window favorable for experimental observation of 
proton emission. 

• Relative energies of different (quasi)particle states 
in medium and heavy mass nuclei are only weakly 
affected by time-odd mean fields. This is because 
additional bindings due to NM show little depen- 
dence on the single-particle state. As a result, the 
present investigation suggests that time-odd mean 
fields can be neglected in the fits of covariant den- 
sity functionals aimed at accurate description of the 
energies of single-particle states. 

• The phenomenon of signature separation [6l| and 
its microscopic mechanism have been investigated 
in detail. It was shown that this phenomenon is ac- 
tive also in two-particle configurations of odd-odd 
nuclei. It is enhanced for configurations having the 
same structure of proton and neutron states; this 
takes place either at ground state or at low excita- 
tion energy in the nuclei at or close to the N = Z 
line. Some configurations away from the N = Z 
line also show this effect but signature separation 
is appreciably smaller. 

The present investigation has been focused on the 
study of time-odd mean fields in the CDFT with non- 
linear parametrizations of the Lagrangian. Point cou- 
pling [65j and density dependent meson-nucleon coupling 
[64J models are other classes of the CDF theories. It 
is important to compare them in order to make signifi- 
cant progress toward a better understanding of time-odd 
mean fields. The work in this direction is in progress, and 
the results will be presented in a forthcoming manuscript 

& 
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